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Model Sampling Applied to the Normal Shock Problem 

by Morris Perlmutter 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Model sampling or monte carlo is applied to the normal 
shock problem. Temperature and density profiles through 
the shock for Mach number of 10 are given. It was assumed 
that the' intermolecular collisions were Maxwellian and that 
the target molecules had a two sided half maxwellian dis- 
tribution. The results show a significant difference from 
previous analytical solutions. Importance sampling is 
used to reduce the solution variance. 


Model sampling analysis has been used previously In 
rarefied gas problems (refs. 1 and 2). This approach is 
extended to the solution of the normal shock problem at 
high Mach numbers for which at the present time no exact 
solution exists. The analytical model is shown in Fig. 1. 
The average quantity Q per molecule transported across 
the scoring cross section p can be obtained by scoring 
the weights times the quantity Q of the sample molecules 
passing that cross section. 



wliere f is the local velocity distribution at point p, 
sp is the number of sample molecules passing p In the 
positive vj direction and W is the weight of the sample 
molecule. The local distribution from which the sampling 
is done ? (v|f/W)^ is obtained by following probable his- 
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tories of sample molecules of weight W from their known 
initial positions to the scoring position. If the weight 
of the entering upstream molecule, Wq, is taken as vp then 
the average Q at position p can be found to be 




'/■pi \ 

[T) m - £*), 


where 


P+O ~ p fo / . . 

- / v 1 >0,x=0 
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% being the maxwellian distribution. The weight of the 
entering upstream sample molecule was chosen so as to mini- 
mize the variance of the density at the upstream scoring 
position. Thus from (Eq. (2)) we can find the local den- 
sity from 



(3) 


Similarly the other local flow characteristics can be 
found. 

The sample histories originate upstream and downstream 
from the shock. The velocity of the sample molecule en- 
tering the first zone upstream from the shock is picked 
from 


f Y l% \ 

V w Jo, ?J>0 
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After the sample molecule enters the zone we must determine 
the path length to collision. The distribution of path 
lengths to collision is given by exp(-&/A )/Ap where the 
mean free path Ap is given by (Kn)p (v s /p). ^ 8 

picked is compared to the distance across the zone, l/p^, 
to determine if a collision has occured. To determine 
the distance across the zones we use 
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If a collision has oc cured we pick a target molecule col- 
lision partner from the distribution of target molecules 
» which was assumed to be a two sided maxwellian as in Ref. 1. 

The parameters of the target molecule distribution axe 
obtained from the sample molecule scoring results and the 
problem iterated. The sample molecule and target molecule 
are assumed to have maxwellian collisions. 

. Results 

The resulting profiles for upstream Mach number of 
10 are shown in Fig. 2. The present results are compared 
to the bimodel model (Ref. 3) and the Navier Stokes results 
given in Ref. 4. The density curve fall along the bimodel 
results upstream of the shock and close to the N.S. down-' 
stream while the temperature curve fall above the bimodel 
and N.S. results upstream and close to the bimodel through 
the rest of the shock. 
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